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Abstract 

'66 , 

, Two-tier networks, comprising a conventional cellular network overlaid with shorter range hotspots (e.g. 

, femtocells, distributed antennas, or wired relays), offer an economically viable way to improve cellular system 

^ ■ capacity. The capacity-limiting factor in such networks is interference. The cross-tier interference between 

o ■ 

macroceUs and femtoceUs can suffocate the capacity due to the near-far problem, so in practice hotspots should use 



(N 



a different frequency channel than the potentially nearby high-power macrocell users. Centralized or coordinated 
frequency planning, which is difficult and inefficient even in conventional cellular networks, is all but impossible 
in a two-tier network. This paper proposes and analyzes an optimum decentralized spectrum allocation policy for 
two-tier networks that employ frequency division multiple access (including OFDMA). The proposed allocation 
Q I IS optimal in terms of Area Spectral Efficiency (ASE), and is subjected to a sensible Quality of Service (QoS) 

requirement, which guarantees that both macrocell and femtocell users attain at least a prescribed data rate. 
^ i Results show the dependence of this allocation on the QoS requirement, hotspot density and the co-channel 
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interference from the macrocell and surrounding femtoceUs. Design interpretations of this result are provided. 



\Q ' I- Introduction 

o ■ 

00 

I Wireless operators are in the process of augmenting the macrocellular network with supplemental 



infrastructure such as microcells [l]-[3], distributed antennas [4] and relays. An alternative with lower 
upfront costs is to improve indoor coverage and capacity using the concept of end-consumer installed 
femtocells or home base stations [5]. A femtocell serves as a low power, short range data access point 
that provides high quality in-building coverage to home users, while backhauling their traffic over the 
IP network. The femtocell radio range (10 — 50 meters) is much smaller than the macrocell radius 
(300 — 2000 meters) [2]. Users transmitting to femtocells experience superior indoor signal reception 
and lower their transmit power, consequently prolonging battery life. The implication is that femtocells 
provide higher spatial reuse and cause less cochannel interference (CCI) to other users. The principal 
arguments in favor of femtocells are summarized in [6]. 
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In a shared spectrum two-tier network (universal frequency reuse), recent research [7] has shown 
that near-far effects arising from cross-tier interference create a fundamental trade-off: Assuming all 
femtocell users seek a common SIR, the product of all Pareto-optimal macrocell and femtocell SIR 
targets is a constant dependent only on the intra- and cross-tier channel powers, irrespective of the 
power control strategy. Perron-Frobenius theory [8] shows that any combination of per-tier target SIRs 
whose product violates the above constraint causes the spectral radius of the normalized channel power 
matrix to exceed unity; hence, there may not be a feasible power allocation for all users. The worst-case 
scenario arise either when a high powered macrocell user on the cell edge causes interference to nearby 
femtocells, or when cell interior femtocell users cause unacceptable interference to the macrocell BS. 
The papers [7], [9]-[ll] have suggested addressing this near-far problem by either 1) forcing femtocell 
users to decrease their target data rates, or 2) employing interference avoidance (e.g. randomized time 
hopping) among both macrocell and femtocell users in order to "avoid" transmitting in the same interval. 
Both approaches have drawbacks because they require either reducing SIR targets at femtocells, or 
accommodating femtocells by altering the transmission scheme at an existing macrocell in order to 
ensure tolerable performance. The subject of this paper is to motivate and propose a spectrum partitioning 
strategy in a two-tier deployment. 

A. The Return of FDMA 

Frequency division multiple access (FDMA)'s resurgence in emerging OFDMA wireless standards 
such as the 3GPP's LTE, WiMAX and 3GPP2's UMB enable the macrocell to perform flexible rate 
assignment [12] across frequency sub bands to users and provide interference management by fractional 
frequency reuse. In femtocell deployments, due to reasons of scalability, security and limited availability 
of backhaul bandwidth, it is reasonable to assume the absence of coordination between femtocells and the 
central macrocell. Further, femtocells are placed opportunistically or randomly by end users. Therefore, 
conventional frequency planning strategies will be very difficult in a two-tier network. 

Assigning orthogonal spectrum resources between the central macrocell and femtocell BSs eliminates 
cross-tier interference. This motivates the orthogonal access spectrum allocation strategy proposed in 
this paper. Next, to avoid persistent collisions with neighboring femtocells in their allotted spectrum, this 
paper proposes that each femtocell accesses a random subset of the candidate frequency subchannels, 
wherein each subchannel is accessed with equal probability. We term this spectrum access strategy as 
F- ALOHA (Frequency ALOHA)^ We motivate F- ALOHA for three reasons. First, F- ALOHA avoids 
transmission delays and increased RF sensitivity requirements for sensing frequency subchannels in the 

'Slotted ALOHA, by convention, implies that the slots are in the time domain 
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presence of channel fading [10]. Next, F- ALOHA provides randomized interference avoidance, since 
neighboring femtocells are unlikely to consistently access identical frequency subchannels. Finally, 
such a transmission strategy offers a decentralized spectrum access by femtocells-eliminating backhaul 
communication cost between femtocells over the intemet-and low complexity. 

F- ALOHA Spectrum Access: If a femtocell transmits over all its allotted subchannels, it may cause 
excessive interference to surrounding femtocells; conversely, accessing only a few subchannels could 
result in a poor spatial reuse. With F-ALOHA, there should be an optimal fraction of spectrum access 
for each femtocell in order to maximize the spatial reuse of spectrum, or in effect the net number 
of simultaneous transmissions per unit area [13], [14]. The spatial reuse is readily expressible using 
the Area Spectral Efficiency (ASE) in b/s/Hz/m^ [15], which is defined as the network- wide spatially 
averaged throughput per frequency subchannel divided by the product of the subchannel bandwidth and 
the area over which the transmissions take place. Based on the stated reasons, we assume (1) downlink 
transmissions from the macrocell and femtocells are frequency orthogonal and (2) femtocells transmit 
using F-ALOHA, and pose the following questions: 

• What is the expected subchannel throughput inside the macrocell [resp. femtocell], as a function of 
interference from neighboring macrocells [resp. femtocells], and terrestrial propagation parameters 
such as path-loss exponent and lognormal shadowing? 

• Given the expected subchannel throughput, the average number of femtocells per cell-site and the 
number of users associated with each BS, how should the bandwidth be partitioned between tiers in 
order to satisfy a Quality of Service (QoS) requirement in each tier? 

• With this spectrum allocation strategy, how much improvement in network-wide ASE does oppor- 
tunistic channel aware macrocell scheduling offer relative to channel blind scheduling? 

B. Related Work 

Existing research on hierarchical cellular systems has mainly focused on channel assignment in two- 
tier macrocell/microcell radio systems [16]-[18]. In the context of this paper, a microcell has a much 
larger radio range (100-500 m) than a femtocell, and generally implies centralized deployment, i.e. by 
the service-provider. These works typically either assume only a single microcell, or regularly spaced 
multiple microcells. Such an assumption may be impractical in femtocell networks because of the 
variations in their placement from one cell site to the next. 

We would also like to clarify the differing objectives between a femtocell and microcell overlay. A 
microcell overlay allows the operator to handoff and load balance users between each tier in response 
to changing traffic conditions [6]. So, the operator can preferentially assign users (e.g high data rate 
users) by converting microcells into data access points [16], [17]. In contrast, femtocells are consumer 
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installed and the traffic requirements at femtocells are user determined without any operator influence. 
In heterogeneous networks, [19], [20] suggest providing vertical handoffs (admission control) and QoS- 
aware routing to switch traffic between cellular and ad hoc technologies. However, these approaches 
may not be viable in two-tier femtocell networks because of the lack of coordination between macrocell 
and femtocell BSs; instead, decentralized strategies for interference management and load balancing 
may be preferred. 

The problem considered in this paper is related to Yeung and Nanda [18], who propose frequency 
partitioning in a microcell/macrocell system based on mobile speeds and the loading of users in each 
cell. Similar dynamic channel allocation schemes have been proposed in [21] and [22]. Their frequency 
partitioning is derived based on choosing handoff velocity thresholds and maximizing the overall system 
capacity, subject to per-tier blocking probability constraints that ignore co-channel interference (CCI). 
In contrast, our work determines the spectrum allocation which maximizes the system-wide ASE 
considering interference from neighboring BSs, path loss and prevailing channel conditions. 

In decentralized networks, Jindal et al. [23] have derived the optimal number of frequency subchannels 
F for a frequency-hopped ad hoc network assuming a fixed data rate requirement per transmitter-receiver 
pair and a target outage probability. Our work, in contrast, assumes a fixed F (as is the case with 
OFDMA), allowing for multiple subchannels to be accessed by each BS, with variable rate transmission 
per accessed subchannel. In a hybrid network composed of ad-hoc nodes and BS infrastructure nodes, 
Zemlianov and de Veciana [24] and Liu et al. [25] have derived asymptotic scaling laws relating to 
the per user throughput as the number of infrastructure nodes increase. In addition to the hierarchical 
nature of our model, the main difference is that this paper assumes single-hop communication. In WLAN 
networks, Bahl et al. [26] have proposed variable center frequencies and variable channel width per 
access point for improving the spectrum utilization and fairness for heavily loaded access points (APs). 
Finally, game-theoretic approaches have been recently investigated for both non-cooperative [27], [28] 
and cooperative [29] spectrum sharing in decentralized networks. 

C. Contributions 

This paper employs a stochastic geometry framework for modeling the random spatial distribution 
of femtocells. Hotspot locations are likely to vary from one cell site to another, and be opportunistic 
rather than planned, so an analysis that embraces instead of neglecting randomness should provide more 
accurate results and more plausible insights. Towards this end, femtocells are assumed to be scattered 
according to a Spatial Poisson Point Process (SPPP) [30], and inter-femtocell interference is modeled as 
Poisson Shot-noise [31]. This model has been used extensively in prior work [32]-[34], and its validity 
has been confirmed in empirical studies [35]. 
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The proposed spectrum allocation maximizes the network-wide ASE in a two-tier network, assuming 
each macrocell transmits to a single user per frequency subchannel, while femtocells access spectrum 
using F-ALOHA. The allocation is determined in two steps. First, the per-tier ASEs are quantified 
based on the propagation environment and neighboring cell interference. With an increasing number of 
macrocell users, the ASE of the macrocell is either fixed (for a channel blind scheduler) or increasing 
(by opportunistic scheduling). For the femtocell network, the ASEs are derived based on the optimal 
spectrum access using F-ALOHA. Next, the optimal allocation is determined as one that maximizes 
the weighted mean of the per-tier ASEs-the weights are given by the fraction of spectrum accessed by 
each tier. The three contributions of this paper are as follows. 

Expected per-tier throughput. The expected per-tier throughput is derived for the macrocell and 
femtocell users accounting for interference from neighboring cells. The maximum ASE of the femtocell 
network is shown to be unchanged with addition of hotspots beyond a threshold. At low femtocell 
densities, the highest femtocell ASEs are attained when each femtocell accesses the entire available 
spectrum. In higher densities, femtocells should use a decreasing fraction of the spectrum; e.g. with an 
average of 100 femtocells in each cell site, each femtocell should access 30% of the available spectrum. 

Spectrum allocation with Quality of Service (QoS). The proposed spectrum allocation maximizes 
the spatial reuse in a two-tier network, subject to a network- wide QoS requirement, which guarantees a 
minimum expected throughput per-user. Differing QoS constraints produce markedly different spectrum 
allocations due to the competing spatial coverage scales in each tier. Notably, a QoS requiring equal per- 
user throughputs in each tier means assigning greater than 90% of spectrum to the macrocell. Conversely, 
an even division of spectrum occurs when the QoS constraints favor femtocells to provide significantly 
higher data rates. 

Scheduling and Spectrum Requirements. Gains generated by channel aware macrocell scheduling 
permeate to femtocells, resulting in a significant spectrum reduction with the proposed allocation. With 
an average of 50 femtocells/cell site and target per-tier data rates of 0.1 Mbps/macrocell user and 10 
Mbps/hotspot user, a channel aware macrocell scheduler provides nearly 50% reduction in necessary 
spectrum compared to a channel blind scheduler. Finally, with increasing number of hotspot users, the 
spectrum requirements in a two-tier network show two extremes. One is a low interference scenario 
where addition of hotspots provides increased spatial reuse, ensuring that the necessary spectrum is 
unchanged up to 110 femtocells/cell site. In a high interference scenario however, the ensuing co- 
channel interference may necessitate a linear increase in required spectrum with hotspot density. 
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II. System Model 

The setup consists of a hexagonal region 7i of radius Rc with a central macrocell BS C providing 
coverage area 1"^^! = ^-R^' which is surrounded by two rings of interfering macrocells. The macro- 
cellular network is overlaid with femtocell hotspots of radius Rf, which are randomly distributed on 
according to a homogeneous SPPP with intensity A/ [30]. The mean number of femtocells per 
cell site is readily obtained as Nf = A/|7Y|. Macrocell users are assumed to be uniformly distributed 
inside each cell site. Femtocells are assumed to provide "closed access" to licensed indoor users who 
fall within the radio range Rf of their respective home BSs^. Let U = Uc + NjUj denote the average 
number of users in each cell site. These U users are distributed into Uc uniformly distributed tier 1 
mobile outdoor users and U f users per femtocell hotspot. 

A. Per-Tier Spectrum Access 

The available spectrum comprises F frequency subchannels each with bandwidth W Hz. We wish 
to determine the optimal partitioning (Fc, Ff), where Fc subchannels are available for macrocell trans- 
missions and Ff — F — Fc subchannels are available for femtocell transmissions. Denote p — Fd F as 
the fraction of spectrum assigned to the macrocell BS with the following key assumptions: 

AS 1: Each femtocell schedules its users in a round-robin (RR) fashion. The macrocell schedules its 
users according to either a channel blind RR or a channel aware proportional fair (PF) scheduler. 

AS 2: The fraction p takes a continuum of values in [0, 1]. For simplicity, the paper does not quantize 
p for having an integer number of frequency subchannels. Consequently, Fc = [pF] and Ff = F — Fc, 
where we use [x] to denote the integer part of a number x. 

Assumption 1 makes it clear that the long term expected throughput per indoor user equals 1/Ufth 
of the average femtocell throughput. The long term expected macrocell throughput is assumed to be 
equally divided among the Uc outdoor users with RR and PF scheduling at the macrocell. With a PF 
scheduler, this assumption is reasonable considering mobility, which ensures that all users receive an 
identical average Signal to Interference Ratio (SIR) over the long term. 

If each femtocell transmits over exactly k frequency subchannels among their allotted Ff subchannels, 
the net portion of accessed spectrum per femtocell equals p/(l — p) where pf = k/Ff. If femtocells 
choose their frequency subchannels independently and with equal probability, F-ALOHA effectively 
"thins" the mean number of interfering femtocells in each frequency subchannel. The probability p of 

^The closed access scheme permits unlicensed users within the femtocell radio range to communicate with the macrocell. 
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a femtocell selecting a given frequency subchannel for transmission is given as: 

Consequently, the set of interfering femtocells per frequency subchannel is a marked SPPP [30] Aj with 
intensity A/p/ = Xfk/Fj. When pf = 1, all femtocells in 0/ access the entire spectrum but mutually 
interfere in all subchannels. For <C 1, femtocells transmit in a small region of spectrum and avoid 
causing mutual interference. This strategy provides a higher spectral efficiency over each frequency 
subchannel, but incurs reduced spectrum utilization because femtocells do not transmit over the entire 
available spectrum. 

B. Channel Model and Variable Rate Transmission 

The downlink channel between each BS and its users is composed of a fixed distance dependent 
path loss, a slowly varying component modeled by lognormal shadowing and Rayleigh fast fading with 
unit average power. For simplicity, thermal noise is neglected at the receiver since cellular systems, by 
nature, are interference limited. 

AS 3: Each user is assumed to track their SIR in each subchannel and feedback the instantaneous 
rate to their BS with zero delay. Further, the channel can support the requested rate as determined by the 
scheduled user with probability 1. Although we acknowledge that imperfect feedback and/or channel 
estimation has a potentially big impact on system capacity, this paper does not account for these effects 
for sake of analytical tractability. 

AS 4: BSs assign equal transmission powers to all subchannels. 
Each BS assigns rate adaptively based on the received SIR per user. Let G denote the Shannon Gap with 
variable rate M-QAM transmission [36]. Assume an instantaneous transmission rate of hi bps/Hz if the 
instantaneous SIR lies in [Fj, Fj+i). Using adaptive modulation with L discrete rates, the instantaneous 
rate Wh in a IV Hz wide subchannel is chosen as: 

h = hi, when SIR e [Fj, F^+i), l<i<L (2) 

hi = log2 [l + ^ bps/Hz (3) 

Assuming identical statistics over all frequency subchannels, the long term expected throughput (in 
b/s/Hz) per macrocell/femtocell in each subchannel is given as: 

L-l 

T^^l- Pr[F; < SIR < Fi+i] + L ■ Pr[SIR > F^] (4) 
1=1 

The expected throughput provided by each macrocell [resp. femtocell] is obtained multiplying the 
expected throughput in (4) by their respective spectrum allocation p [resp. — p)]. 
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III. Spectrum Allocation and Per-Tier Expected Throughputs 

Let spectrum WF be partitioned such that the macrocell BS transmits over a portion p, while femtocell 
BSs transmit over the remaining 1—p fraction of the spectrum. Let Tc{p, Uc) be the long term throughput 
(in b/s/Hz) in each subchannel provided by the macrocell^. Obtaining Tc requires calculating the average 
rate per subchannel in (4) after spatially averaging the SIR over all locations, and accounting for the 
interference from two rings of transmitting macrocells. 

Let each femtocell access a portion p/ of its allotted spectrum using F- ALOHA, servicing its users 
in a RR schedule. Define T/(p/A/) as the expected femtocell throughput in each frequency subchannel, 
which is determined by the intensity pf\f of the marked SPPP Aj. With universal frequency reuse 
across all macrocells, the area spectral efficiency (ASE) of the macrocell (resp. femtocell) network is 
given as: 

\rL\ \ri\ 
The factor NfPf represents the mean number of transmitting femtocells in each subchannel. With 
bandwidth W, the per-tier throughputs (in b/s) per subchannel can be calculated by multiplying the 
ASEs in (5) by The network- wide ASE is therefore given as: 

ASE = pASEc + (1 - p)ASE/ 

= ^[pTe(p, t/e) + (1 - p)NfPfTf{pf\f)] (6) 

The expected network throughput (in b/s) over the WF wide spectrum is obtained by multiplying (6) 
by WF\H\. Before determining the spectrum allocation, we first stipulate a QoS requirement 77, which 
ensures that users in either tier are guaranteed a minimum expected throughput. By implication, rj also 
regulates the maximum amount of spectrum that any tier can receive. 

Definition 1: The QoS parameter 77 guarantees that the expected throughput per user in one tier is at 
least 77/(1 — 77) w.r.t the other tier. 

Choosing different 77 enables assigning different priorities (QoS) to one tier relative to the other. For 
example, setting 77 = 0.5 ensures that users in both tiers obtain identical expected rates. On the other 

hand, decreasing 77 favors assigning greater spectrum to the tier providing a higher expected throughput 
per active user. 

Given a total available spectrum of 1 Hz, the problem is to determine the optimal spectrum allocation 
p over all possible spectrum partitioning strategies u e [0, 1] between the macrocell and femtocells. 



'The use of Uc and p within parenthesis is to account for a macrocell scheduler which can provide diversity gains by scheduling 
users according to their channel variations [37] 
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The proposed spectrum allocation maximizes the network-wide ASE with a QoS constraint 77 on the 
minimum expected per-tier throughput/user, as shown below: 

P=\h\ ^n{u;, U,) + (1 - u;)NfPfTf{pfXf) (7) 

subjectto mm{r,,„(ij),r^,„(a>)}>T;(r„(w)+r,.„(w)) (8) 

where T„(.) 4 ^^I^ and r,„(.) 4 d 

Here Tc„(a;) and Tf^u{u!) are the expected throughputs for a macrocell and femtocell user respectively. 
Whenever the average subchannel throughput Tc{oj, is independent of the spectrum uj assigned to the 
macrocell, the objective function in (7) is an affine function w.r.t u. The following proposition derives 
the optimizing p considering that the maximum is attained at the extremal points of the constraint set: 
Proposition 1: If the expected macrocell throughput per subchannel is independent of the total 
spectrum allocated to the macrocell u, i.e. Tc{uj,Uc) = Tc{Uc) "iuj G [0,1], the optimizing p in (7) 
satisfies the QoS constraint with equality, belonging to a set with two candidate spectrum allocation 
assignments given as: 

p* e {x, 1 — x}, X = 

Proof: Since Tc{u, — Tc{Uc) Vo; e [0, 1], the optimization problem in (7) is to determine the 
optimal p which maximizes a convex combination of Tc and NfPfTf{pfXf) with a linear constraint. 
Consequently, the argument maximizer is located at the extremal points of the constraint set (8). Solving 
for the p which satisfies (8) with equality yields (9). ■ 
Remark 1: Without a QoS requirement (allowing r] 0), the objective function in (7) is a convex 
combination of the macrocell and femtocell throughputs which is maximized at the extreme points 
p e {0, 1}. Such a partitioning is clearly unfair since it results in a greedy allocation of the entire 
spectrum to one tier. 

For a generic macrocell scheduler-when Proposition 1 may be inapplicable- (7) is a one dimensional 
optimization problem that can be solved efficiently for a given 77 using a numerical search. 

A. Macrocell Throughput: RR Scheduling 

Assuming that the central macrocell BS Co in the hexagonal region H is placed at the origin, the 
normalized positions of the interfering BSs e B, /c = 1 . . . 18 are represented in polar form bk, k E B 
using MATLAB notation as: 



1 + 



V Uc PfTfipfXj 



(9) 
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Let ho denote the Rayleigh fast fade (with exponentially distributed channel power |/ioP with unit mean) 
between the central macrocell BS Co and its scheduled user 0. Denoting the Euclidean norm by ||H, 
the expression for the received SIR for macrocell user at position r is given as: 

eo 



SIRc(r) 



where ^/(r) = 

keB 



\ho\ \\r/Rc 
r - bfc 



Rr 



OA; 



(11) 

(12) 



Here etc represents the outdoor path loss exponent and |/iofeP ~ exp(l) is the exponentially distributed 
channel power between interfering BS and the user of interest. The RV ©o [resp. ©ofe] is the 
lognormal shadowing between the central BS [resp. interfering BSs] and the desired user, which are 
distributed as LN(CyU.c,dB, C^crc,dB)' where ( = 0.1 In 10 is a scaling constant. For analytical tractability 
in the reminder of the paper, our paper makes the following assumption regarding the distribution of a 
composite lognormal-exponential RV: 

AS 5: The distribution of a composite lognormal-exponential RV ©fe|/iofeP is modeled as a lognormal 
distribution using Turkmani's approximation [38]: 



1 



xctjV 27r 



exp 



- (In a; — /Xj 



(13) 



f^i = C(/^c,dB - 2.5), (J, = C7<Tc'dB + 5.572 (14) 
AS 6: For a fixed r, using the moment generating function based technique described in Mehta et 
al. [39], the sum of |B| independent, but not identically distributed lognormal RVs in the expression 
^ = Z^fcei^dl^ ~ /^c)"""^*: is approximated by a single lognormal RV with parameters 

LN(;.,(r),^Kr)). 

Using Assumptions 5-6 and (11), SIRc(r) is distributed according to a lognormal RV LN(/ic(r), (T^(r) 
where iic{^) — iJ'si'i^) — l^ii'"^) and ac{r) — ^/ag(rY+~aJ(r). Then, the distribution of the SIR for a 
mobile tier 1 user at position r w.r.t the central BS is given by 



Pr[SIRc(r) <r|r] = 1-Q 



ln(r||r/i?,|r=)-/ic 



(15) 



where Q,{x) = JT' ^ ^^^^dt is the complementary cumulative distribution function (CCDF) of a 
standard normal. Defining a(r) = ^"^"^^^'■'"^ and b = -^^pj, (15) simplifies to 

Pr [SIRc(r) < r|r] = 1 - Q[a(r) + 6(r) In | \r/R,\\] (16) 

Averaging (16) over a hexagonal cell region is difficult. Alternatively, the spatially averaged CDF of 



SIRc can be obtained approximately by considering an circular region of radius 



n 



, which results 



in the same area as the cell site H. To calculate the spatial throughput inside this circular region, the 



11 



paper divides the region into M non-overlapping annuli. For tractability, a simplifying assumption is 
that all users inside an annulus experience identical shadowing statistics (i.e. identical /xc(r) and ac{r)). 
Denoting the distance of the user from Co by | |r| | = i?, the following lemma derives the expected spatial 
throughput by averaging SIRc(-R) inside a circular annulus with inner radius -Ri and outer radius R2. 

Lemma 1: The spatially averaged SIR distribution inside a circular annulus with inner radius Ri 
and outer radius R2 is given as: 

ErIFv (SIRe < T\R, <R< R2)] = 1 - T52^-B2^[RlC{a2, h) - RlC{a,, b)] (17) 
Where C{a, h) ^ Q(a) + exp (|^^ j Q (|^^ j (18) 

^ = r5~\ — — r5~\ (1^) 

ac[R2) crc[R2) 
02 = a + 6 In {R2/Rc) ,ai = a + b\n (Ri/Rc) , (20) 
Proof: See Appendix I. ■ 
Lemma 1 provides a simple method for estimating the cell-averaged macrocell throughput per sub- 
channel. The probability that a user lies in an annulus with inner radius Rm-i and outer radius R^ 
(1 <m < M with Ro — 0) equals — . We make use of assumptions 5 through 6 for computing 

the shadowing parameters ac and fic at discrete locations /?„i, 1 < m < M where Rm = \f^-- The 
spatially averaged SIR distribution for a macrocell user is therefore approximated as follows: 

Pr (SIRc < r) = Efl[Pr (SIRe(i?) < T)] (21) 
E«[Pr (SIRe < r|i?^_i <R< R^)] ■ ^ ""^^ (22) 



= 1 - TTTr*^ (ai + bi In — bi 



1^1 V R< 



Rl^C (am + bmin^,bmj -i?^_iC fa„ + 6^1n 6^+1 



m=2 



\H\ 



''C 



(23) 



where = 7^-^^ and b^ = ac/ac{Rm) 

where (22) approximates (21) by spatially averaging SIRc over M different annulus. Equation (23) is 
obtained by substituting (17) inside the conditional expectation in (22) and the corresponding probability 
that the user lies in annulus m,! <m < M. Combining equations (4) and (21), the average macrocell 
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throughput Tc in a given subchannel is expressed as: 



L-l 



< siRc(i?) < r^i)] + L ■ E^[Pr(siRe(i?) > r^)] 



1=1 

L-l 

= ^ / ■ (E«[Pr(SIR,(i?) < r^+i)] - EK[Pr(SIR,(i?) < Ti)]) + L ■ ER[Pr(SIR,(7?) > T^)] (24) 
1=1 

Figure 1(a) plots Tc (in b/s/Hz) with RR scheduling as a function of the outdoor path-loss exponent etc 
for the system parameters in Table I. The close agreement between theory and numerical simulations 
indicates that the theoretically obtained SIR distribution is an accurate approximation for practical 
throughput in a macrocellular environment. 

B. Macrocell Throughput: PF Scheduling 

In contrast to a RR scheduler, a PF scheduler enables macrocell users to compete for resources based 

on their requested rates normalized by their average throughput thus far. Consequently, the macrocell 

selects the user with the highest rate relative to their average rate. During the transmission interval n in 

subchannel m, denote Rk[m, n] as the requested rate for user /c, 1 < A; < Uc, located at position w.r.t 

the central macrocell C. Let Rk\n] be the windowed mean throughput obtained by user k over the Fc 

frequency subchannels allocated for macrocell transmission. The PF scheduler selects the user k whose 

current supportable rate is highest relative to their mean rate. The scheduling policy per subchannel m 

with equal per-subchannel transmission powers (Assumption 4) is described as: 

l( \ Rk[m,n] 
l<fe<?7c Rk[n\ 

Note that mobile user k calculates Rk[m,n\ using (2) and (11) respectively. The windowed throughput 
per user prior to transmission interval {n + 1) is updated according to the rule: 

1-1 

Rk[n + 1] = (1 - -)Rk[n] + ^ J] ^^C"'' ^l^C^ = ^)\^<k<Uc (26) 

m=l 

where ![■] is the indicator function determining whether user k is scheduled during transmission interval 
n in frequency link or not. The window size N is di parameter that is selected considering the delay 
tolerance for each user. Choosing a smaller N enables a given user to be scheduled more often, 
whereas choosing larger N relaxes the fairness constraint and allows the scheduler to wait longer 
before scheduling a user. By the strong law of large numbers, the average throughput per frequency 
subchannel for a given set of user positions is obtained from the sample average over a long duration 
and expressed as: 

E[i?(F„t/e)|ri,---rc;J= lim - V V G {1, 2 ■ • • f/e} (27) 
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where the expectation on the left hand side is over the joint pdf of all channel gains between users and 
their serving and interfering BSs. The spatial averaged subchannel macrocell throughput is obtained by 

averaging (27) w.r.t the joint pdf /ri,...Rj^^(-) and given as: 

Te(p, Uc) = Er,,..r^JE[^(Fc, C/e)|Rl = T,, • • • R[/, = Ft/J] (28) 

Using (28) to compute Tc(p, Uc) is analytically intractable. This paper resorts to numerical simulation 
to empirically estimate Tc{p, Uc), which is used to derive the bandwidth partitioning. In the simulation, 
the number of subchannels is set as Fc = 1 with a link bandwidth W — 15 KHz and a PF window 
parameter N — 500 OFDM symbols. Each mobile is moving atv — 13.34 m/s (30 mph) and the per-link 
throughput (27) is averaged over 500 drops, with 8000 trials/drop for modeling time-varying Rayleigh 
fading. The Rayleigh fading is held fixed over a duration Tc = OA/fd where fd = is the doppler 
frequency at a carrier frequency fc = '^ GHz. Figure 1(b) compares the performance of PF (numerically 
evaluated) versus RR scheduling for different Uc- Exploiting channel variations through proportional 
fairness roughly doubles the expected subchannel throughput. 

C. Femtocell Throughput 

Since femtocells are modeled as randomly placed on according to a SPPP $/ with intensity A/, the 
interference experienced by a femtocell user depends on the distances of these interfering BSs relative to 
the user and their respective channel gains. Using F-ALOHA, the interfering femtocells form a marked 
SPPP Af C with intensity PfXf. In a given frequency subchannel, the cochannel interference Ifj 
experienced by a user within femtocell Fq is given as: 

If J = ^/0ofe|/iofc|'kofc|-"^ (29) 
fceA/ 

where ©ofe ~ LN((^yUjo dB, C^o-j^^jg) represents the lognormal shadowing from femtocell to user 
and |/iofeP is the exponentially distributed channel power between interfering femtocell F^ and user 
inside Fq. Denoting the exponentially distributed channel power between Fq and user as |/ioP> the 
received SIR is given as: 



SIR. = 



BfQol 


^0 


|2 


\Rf 


\-l3f 




Ok\ 


hok\ 


2j 







(30) 



fceA^ 



Here, user is assumed to be on the edge of the home femtocell Fq and xok represents the locations of 
the interfering femtocells Fk w.r.t user 0. The term ©o ~ ^^{Cl^fi,dS: C^<^/i,dB) is the indoor lognormal 
shadowing, and — ©ol^oP [resp. ^^ofe — ©ofcl^oitp] are the effective channel gains from the 
desired [resp. interfering BSs]. The terms o;/ [resp. /?/] represent the path-loss exponents resulting from 
interfering transmissions [resp. in-home transmissions] to the user of interest. A simple model is used 
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to distinguish between the fixed losses arising from in-home and interfering transmissions: specifically, 
home users are insulated against interfering femtocell transmissions through double penetration losses 
arising from external wall partitions [40]. Consequently, Af and Bf (in dB) are related as A/,dB = 



-B/,dB + 2P/,dB where Pf 



is the wall penetration loss. 



Using Assumption 5, the channel gain ©ol^oP is well approximated as a lognormal rv 
LN(//5, (7|). Similarly, the channel gains ©ofel^ofepVA; are approximated as iid rv's distributed as 
LN(/X7,(7|). Equation (30) then simplifies to: 



*o|i?/| 



keAf 



(31) 



The closed form distribution of the Poisson SNP Ifj = Yl Pf^ok\xok\ "^^ is known only when af — A 

keAf 

[31]. However, tight lower bounds on Pr {Ifj > y) are obtained by only considering femtocells whose 
interference individually exceeds y. Using this idea, the following lemma, derived in Weber et. al. [14, 

Theorem 3] provides an asymptotically tight lower bound on the tail distribution of Ijj. 

Lemma 2: [14, Theorem 3] With randomized transmissions and lacking power control, the lower 
bound on distribution of Ifj is given as: 



Piilfj >y)>l-exp [-nXfPfE[^'/]Ppy-'f] 



(32) 



where 5f = When af — A, Ifj is distributed as: 



Pr {Ifj >y)^l- erfc I ^^^^^ ^ 



(33) 



Lemma 2 provides the relationship between the density A/p/ of interfering femtocells in A/ and the 
distribution of the CCI at a femtocell. For fixed y, as pf — > 0, the tail probability Pr(//j > t/) — > in 
(32) indicating that selecting fewer subchannels using F- ALOHA transmission provides greater resilience 
against persistent collisions from nearby femtocells. The distribution of the femtocell SIR in (30) is 
obtained as: 



Pr (SIR/ < r) = Pr 



— 



^o\Rf\~^f 



E Pf^ok\xk\-"f 
•■keAf 



< r 



Pr 



Ok\Xk 



\keAf 



> 1 — < exp 



-TTXfPfEi^f 




^l-E^,[exp{-pfKfT'f%'^: 
where, Kf = 7rXfE[^^/]{P]\Rf\^fYf 



(34) 
(35) 

(36) 
(37) 
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where (35) and (36) follow by conditioning on ^^o^ assuming independence of and ^^ofe VA; e Aj, 
and applying (32). Although it is not possible to obtain a closed fonn expression for the expectation 
in (37), the distribution of SIR/ can be calculated numerically. The mean subchannel throughput Tf is 
calculated by combining (4) and (37): 

L-1 

Tf{Pf\f) = ^ ^ Pr(ri < SIR/ < r^+i) + l ■ Pr(siR/ > r^) (38) 

1=1 

L-1 

E^Jexp {-pfKfV\!^,%'^) - exp {-pfKfT'/%'^)] 

1=1 

+ L ■ E^Jexp {-pfKfV'l%'^)] (39) 

The approximation in (38) is because the right-hand side in (37) is a lower bound on Pr(SIR/ < V). 
Figure 2 plots the femtocell throughput (1 — p)pfTf (in b/s/Hz) assuming the entire bandwidth is 
allocated to femtocells (p = 0). Black colored curves plot results of numerical simulations. Two cases 
are considered namely (1) high attenuation (marked "HA" with af — 4:,Pf^dB — 10) and (2) low 
attenuation (marked "LA" with a/ = 3.5, Pf,dB — 2) from neighboring femtocells. Setting Pf — 1 and 
assuming Nf — 50 femtocells/cell site, the femtocell throughput falls from approximately 4.5 b/s/Hz in 
a HA environment to nearly 0.5 b/s/Hz in LA scenario, indicating the sensitivity of femtocell throughput 
to propagation from nearby femtocells. 

To calculate the optimum p/, we resort to maximizing the ASE per subchannel. This is analogous 
to answering the question: What fraction of subchannels should each femtocell access to maximize 
spatial reuse? At this critical p/, the F- ALOHA access by each femtocell is optimally traded off 
against neighboring femtocell interference in each subchannel. Mathematically, p/ is the solution to the 
following optimization problem: 

p} = A/argmaxer/(eA/) (40) 

ASE} = p}A/T/(p}A/) (41) 

To justify (40), observe that there are an average of |7Y|p/A/ transmitting femtocells per subchannel. 
With F- ALOHA access of < ^ < 1, each femtocell obtains an average subchannel throughput of 

Tf{9), resulting in ASE/ equaling \f9 ■ Tf{9Xf). Alternatively, given any allocation p, (40) computes 
the F-ALOHA spectrum access p/ which maximizes the mean overall throughput (1 — p)pfTf{pfXf) 
per femtocell. 

Remark 2 (Boundedness of the ASE): The ASE in (40) depends on the effective intensity Xf9 of 
interfering femtocells per subchannel. With increasing A/, provided p/ < 1, then the intensity of A/ 
given as A/p/ remains constant, implying that the optimal p/ is a monotone decreasing function of A/. 



16 



Consequently, if p/ < 1 for a given A/, the maximum ASE per subchannel is fixed. This also means 
that with increasing A/, the network-wide femtocell throughput equaling |?i|VFF • (1 — p)ASEj grows 
linearly with (1 — p). 

Fig. 3 plots (40) for different Nf with o;/ = 3.5 and -P/,ds = 2. In all cases, the highest ASE is fixed 
at nearly 0.000121 b/s/Hz/m^ vaUdating Remark 2. With a low femtocell density {Nf = 10), the best 
strategy is to access the entire spectrum from bandwidth partitioning. In a dense network {Nf — 100), the 
ASE is maximized when each femtocell accesses approximately 30% of the available spectrum. Further, 
in (40), as long as = 1, each femtocell accesses the entire available spectrum (1 — p), consequently 
Tf decreases with addition of femtocells. However, if p*j < 1, Tf = ASE//(A/pj) remains constant 
with increasing Xj (Fig. 3(b)). However, as Aj ^ oo, since p e (0, 1), the mean overall throughput per 
femtocell approaches zero, as the following limit shows: 

lim Tf{l - p)pf < lim TfPf = (42) 

A/— »oo Ay-+oo 

One may explore the dependence of the mean overall femtocell throughput Tfpf{l — p) on the spectrum 
allocation p and F-ALOHA access pf. Equivalently: With increasing femtocell density A/, can increasing 
allocated spectrum (1 — p) to femtocells counterbalance decreasing F-ALOHA spectrum access Pf to 
result in a higher mean femtocell throughput? 

This question is answered by the following condition: Given an allocation pi at femtocell density 
A/, let Tf^i and pf^i be the mean subchannel throughput and the optimal F-ALOHA access respectively. 
On increasing A/ by 5\f with allocation p^, let the corresponding quantities equal T/,^ and pf^h- The 
femtocell network is defined diS fully-utilized [ resp. sub-utilized] if a marginal increment in the femtocell 
density reduces [resp. improves] the mean throughput per femtocell as given below: 

(1 - Pi)pf,iTf,i ^ (1 - Ph)pf,hTf,h 



Xf Xf + SXf 

^ ^~ < ^^^f'l . + = IjlPf£ (43) 
l-p/^ASE/,/j Xf Tf^hPf.h 

Equation (43) reflects the competing effects of increasing allocation (1 — p) and decreasing F-ALOHA 

access p/ (or increasing Xf) in determining the net femtocell throughput. 

IV. Numerical Results 

Results are presented in Figs. 4 through 7 with the system parameters in Table I. The number 
of users in each tier is controlled by varying Nf. To model varying data-rate requirements inside 
femtocells relative to the central macrocell, QoS values of 77 = 0.5 (equal per-user throughputs in each 
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tier) and rj — 0.01 (favoring lOOx higher throughput/femtocell user relative to macrocell users) are 
considered. Two propagation scenarios are presented namely 1) High Attenuation (HA) of neighboring 
femtocell transmissions-modeling low interference-with parameters a/ = 4 and Pf,dB = 10 and 2) Low 
Attenuation (LA) scenario by setting aj = 3.5 and P/^b = 2. 

Fig. 4 shows the allocation using (9) with RR scheduling per-tier. With equal average per-user 
throughputs (rj — 0.5), nearly 90% of the overall bandwidth is assigned to the macrocell. The central 
macrocell serves a higher number of users, who experience relatively poor reception. Equalizing per-user 
throughputs consequently requires a significantly higher allocation to the macrocell. As rj decreases, 
femtocells require more spectrum for providing greater indoor capacity; eg. in a LA scenario with 
7] = 0.01 and Nf = 50 femtocells/cell site, nearly 70% of spectrum is allocated to femtocells. 

Fig. 5 plots the ASEs of the two-tier network using (6). In a LA scenario with rj — 0.01, the high 
degree of co-channel interference results in the ASE maximized with fewer than Nf — 50 femtocells. 
Following Remark 2, this indicates that adding more femtocells does not provide additional spatial reuse. 
In all other cases, the ASEs monotonically increase with Nf indicating increasing spatial reuse with 
addition of femtocells. To show benefits of opportunistic scheduling, a PF scheduler provides nearly 
2.3x [resp. 1.35x] ASE gains relative to a RR scheduler in a HA scenario with QoS parameter rj = 0.5 
[resp. rj — 0.01] and Nf — 110 femtocells/cell site. 

Fig. 6 plots the expected throughput per femtocell {l—p)pfTf as a function of Nf and rj. For 77 = 0.5, 
the throughputs monotonically increase with Nf indicating that increasing spectrum allocation (1 — p) 
counteracts the effects of decreasing p/T); in effect, the femtocell network is sub-utilized. With rj — 0.01 
in a LA environment however, the femtocell throughputs decrease with increasing Nf, indicating that 
the femtocell network is fully-utilized. 

Fig. 7 plots the minimum required spectrum WF, which satisfies a target average data rate of = 0.1 
Mbps for each macrocell user. For each femtocell, corresponding to QoS parameter values 77 = 0.5 and 
rj — 0.01, we consider target average data rates Df — Dc{l — r])/rj equaling 0.1 and 10 Mbps/user. 
Since Proposition 1 ensures that the QoS constraint in (7) is binding, the required spectrum WF is 
given as: 

WF = = ^( Df (44) 

pTc (1 - pjpfTf 

Two key observations are: First, a channel aware macrocell provides significant savings in the spectrum 
necessary to meet Df and D^, eg. with rj = 0.01 and Nf = 50 femtocells/cell site in a HA scenario, the 
spectrum reduction is nearly 50% (10 MHz). Next, spatial reuse and spectrum requirements with the 
addition of femtocells are markedly different depending on attenuation from neighboring femtocells. For 
example, in a LA [resp. HA] scenario, the spectrum requirement WF increases [resp. decreases] with 
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increasing hotspot density indicating the femtocell network is fully-utilized [resp. sub-utilized] with the 
per-tier spectrum allocation in (9). 

V. Conclusions 

This paper has proposed a decentralized spectrum allocation strategy as an alternative to central- 
ized/coordinated frequency assignment in a two-tier network. The proposed allocation depends on the 
per-tier throughputs, the loading of users in each tier and the QoS requirements, accounting for co- 
channel interference and path-losses. With a randomized spectrum access strategy, femtocells should 
access a decreasing fraction of their allocated spectrum with increasing femtocell density, in order to 
maximize spatial reuse. Spatial reuse benefits derived from channel aware macrocell scheduling result in 
nearly 50% spectrum reduction for meeting target per-tier data rates. In a low interference scenario where 
addition of hotspots provides increased spatial reuse, the spectrum requirement is unchanged up to 110 
femtocells/cell site. On the other hand, the limited spatial reuse in high interference scenarios necessitates 
increasing spectrum with addition of femtocells. These insights provide guidelines on performance of 
decentralized spectrum allocation in the two-tier networks. 

Appendix I 

Inside a circular annulus of small width, one can assume that a user experiences identical shadowing 
statistics from interfering BSs. For convenience, the parameters a and b are chosen when the user is at 
the outer edge (R — R2) of the annulus, as shown in equation (19). Given Ri < R < R2 with the pdf 
function fR{r\Ri < R < R2) — ^a^^a , the spatially averaged outage probability is obtained as: 

Eflpr (SIRc < T\Ri <R< R2)] = 

1 - EulQia + 6 In R/R^) \ Ri < R < R2\ (45) 

= 1 - p2 ^ p2 Q[« + ^ In r/R,]rdr (46) 

1 { /"^^ 2 

J, ^Q[«2 + Mnr/i?2]rdr 

- Rl -^QWi + Mnr/i?i]rdr| (47) 

where (47) is obtained by substituting a2 = a + bin R2/ Rc and ai = a + blnRi/Rcin (46). Combining 
the definitions in (18) and (20) with the identity [36, Pg. 55]: 

2 

C{a, b)^—^ Q(a + 6 In r/R)rdr (48) 

R Jo 

and plugging into (47), the result follows. 
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TABLE I 

System Parameters 



Symbol 


Description 


Value 


Rc, Rf 


Macrocell/Femtocell Radius 


288 m, 40 m 


U 


Total users per cell site 


300 


Uf 


Users per femtocell 


2 


Pf,iB 


Wall penetration loss 


2 dB, 10 dB 


G,L 


Shannon Gap, Modulation Levels 


3 dB, 8 




Path-loss exponent (Macrocell Outdoor) 


4 


otf 


Path-loss exponent (Femtocell to Femtocell) 


3.5,4 




Path-loss exponent (Inside Home Femtocell) 


3 


C"c,dB, 0"/i,dB, ""/o.dB 


Lognormal Shadow Parameters 


8 dB, 4 dB, 12 dB 




Path Loss Exponent [a^] Users Per Macrocell 

(a) Tc versus outdoor path-loss exponent Qc (b) Tc with RR and PF Scheduling, ac = 4 

Fig. 1. Spatially averaged macrocell subchannel throughput Tc (b/s/Hz) 
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(b) Femtocell ASEs and subchannel throughputs 



Fig. 3. Femtocell Area Spectral Efficiencies ^^.^^ 
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Fig. 4. Optimal spectrum allocation p for varying QoS parameter r) 
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Fig. 5. Area Spectral Efficiencies in a two-tier network for varying QoS parameter 77 
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Fig. 6. Average network- wide throughput pfTf{l — p) (in b/s/Hz) provided by femtocells in their allocated spectrum (1 — p) 




Fig. 7. Required spectrum WF meeting a target average data rate of Dc = 0.1 Mbps for each macrocell user, given Round-Robin and 
Proportional-Fair scheduling at the macrocell. 



